Exciting a d-density wave in an optical lattice with driven tunneling 
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Quantum phases with unusual symmetries may play a key role for the understanding of solid state 
systems at low temperatures. We propose a realistic scenario, well in reach of present experimental 
techniques, which should permit to produce a stationary quantum state with d x 2_ y 2 -symmetry in a 
two-dimensional bosonic optical square lattice. This state, characterized by alternating rotational 
flux in each plaquette, arises from driven tunneling implemented by a stimulated Raman scattering 
process. We discuss bosons in a square lattice, however, more complex systems involving other 
lattice geometries appear possible. 

PACS numbers: 32.80.Pj, 03.75.Hh, 03.75.Nt, 11. 15. Ha, 75.10.Jm 



Solid state systems at low temperatures may ex- 
hibit exotic quantum phases, which involve order pa- 
rameters with unusual symmetries. Various researchers 
have considered order parameters with d-wave symme- 
try 0) S B HL 111 with the hope of shedding light on 
the yet unexplained pseudogap phase in cuprate super- 
conductors. Quantum states with d 2 .2_ y 2-wave symmetry 
(staggered flux states), characterized by alternating rota- 
tional flux in each plaquette, have been shown to arise in 
the context of the two-dimensional Hubbard-Heisenberg 
model p|, or in the context of the t — J model [2j. Un- 
fortunately, decisive experiments, which could clarify the 
role of such unconventional symmetry broken phases, are 
hampered by the effects of inevitable disorder in solid 
state systems. 

In recent years, atomic or molecular quantum gases Q 
loaded into the perfectly periodic light-shift potentials 
of optical lattices have emerged as synthetic laboratory 
systems which permit to experimentally simulate basic 
models of low temperature physics while keeping all sys- 
tem parameters under control |7fl. A seminal example 
is the observation of the superfluid-insulator phase tran- 
sition within the Bose-Hubbard- model [1, 0, Ever 
since there have been rapidly growing activities to en- 
gineer optical lattice models of more involved cases like 
paired fermionic systems ll|, 1 1 21 ] , quantum Hall systems 
[l3|, or to study unusual lattice geometries like Kagome 
lattices p~3|. According to a recent proposal [Hj], quan- 
tum states with unconventional order could be imple- 
mented in an optical lattice of bosons if their internal 
structure allows for additional engineered interactions. 
Unfortunately, the existing proposals to extend the scope 
of optical lattices to the exciting realms beyond the con- 
ventional Hubbard model are experimentally extremely 
demanding or even involve requirements which cannot be 
met with state of the art technology. 

In this article we discuss a simple realistic scenario 
which allows for the experimental implementation of an 



order parameter with d-wave symmetry in a two dimen- 
sional square optical lattice of simple polarizable bosons. 
Our proposal relies on the existence of a metastable sta- 
tionary solution of the Gross-Pitaevski (GP) equation 
with d x 2_ y 2-w&ve symmetry. In order to excite this 
state, we propose a method based upon Raman scat- 
tering. Apart from a non-zero polarizability, no inter- 
nal structure of the bosons is required, a circumstance 
which should dramatically facilitate experimental inves- 
tigations. We exploit the interference between two ap- 
propriate light fields with slightly different frequencies, 
which produces a two-dimensional square array of micro- 
scopic rotors. These rotors generate angular momentum 
with alternating rotational direction for adjacent plaque- 
ttes. In the framework of the Bose-Hubbard-model, we 
show that this Raman coupling leads to driven ring-like 
tunneling characterized by anisotropic time-dependent 
hopping. This opens up a new dimension of phase space, 
possibly giving rise to new quantum phases, accessible in 
the clean environment of an optical-lattice system. Here, 
we restrict ourselves to bosonic atoms in a square lat- 
tice geometry, however, more involved scenarios are in 
reach of present experiments. For example, a molecular 
optical lattice could be used to access the BCS regime 
[12|. D-density waves subjected to triangular or quasi- 
periodic lattices should provide exciting new options to 
study frustration phenomena Q. 

Let us consider solutions of the time-independent GP- 
equation in an external potential Vt rap (r) 



A + U trap (r)+ ff p(r) V(r) =/x^(r) (1) 



2m 

of the form ip(r) = a(r)<fi(r), with p(r) = IVK 1 *)! 2 the 
particle density and \x the chemical potential. The col- 
lision parameter g = Ai:h 2 a/ra, with m the atomic mass 
and a the s-wave scattering length, is positive for re- 
pulsive collisions. Upon use of the eigenfunction of the 
kinetic energy <p{x,y) = e l6 / 2 svn{kx) + e~ ie / 2 sin(fcy), 
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-{h 2 /2m)A<j) 



E rec (j> with the recoil energy 
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H 2 k 2 /2m, Eq. (1) becomes -(h 2 /2m)[Aa/a 



trap 



99 = 9- 



Assum- 



2Vln(CT)Vln(0)] + E Te 
ing a positive, slowly varying envelope function cr(r) 
such that |Vit| <C fccr and |Acr| <C fc 2 c, one may 
approximate \Aa/a + 2Vln(<r)Vln(0)| < \Aa/a\ + 
2|Vln(ff)||Vln(0)| < |Acr/cr| + 2fc|Vln(cr)| < 3fc 2 and 
thus |Act/ct + 2Vln(cr)Vln(0)| < £; rec . Consequently, 
the expression between square brackets may be neglected 
as compared to E rcc . This leads to E rcc + Vtrap + 99 = 9 
with the possible solution 



fi = E Tl 



Vtrap — 



'99 , 



(2) 



where p = a 2 \4>\ 2 and \4>\ 2 = sm 2 (kx) + sin 2 (ky) + 
2cos(#) sin(fcx) sin(fcy). With = |0|e iS, ^ ,3/ ) the phase 
S(x, y) becomes 



S(x, y) = arctan 



tan 



sm(kx) — sin(fcy) 



+ - sgn[3(0)] (1 



sin(fea;) + sin(fcy) 
sgn[5R(0)]), 



(3) 



where sgn[z] denotes the sign of z. The particle cur- 
rent density j(x,y) = \4>(x, y)\ 2 v(x, y) connected to 
the velocity field v(x, y) = (H/m)'VS(x,y) evaluates to 
j(x,y) = (H/m) sin(0) V x z sin(fcx) sm(ky), i.e., it is 
purely vortical (see Fig. 1(a)). In the case sin(0) = 1, 
where j is maximal, the local particle density scales as 
\4>{x,y)\ 2 — sin 2 (fcy) + sin 2 (/cx) (see Fig. 1(b)) and the 
velocity field satisfies 



h V x z sm(kx) sin(fcy) 

v(x,y) = it, — — o-r 

m sin [ky) + sin (kx) 



Vx w(x,y) = 



(4) 
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According to Eq. (jU, the structure of v(x, y) is that 
of a vortex-anti-vortex lattice with vortex-filaments at 
positions kx,ky G 7rZ. At the potential minima (white 
regions in Fig. 1(b)) the particle density is maximal and 
v(x, y) exhibits quadrupole symmetry. The vortices are 
pinned at the potential maxima of the lattice (black 
regions), in accordance with results obtained for plain 
vortex lattices prepared in large scale traps and subse- 
quently exposed to an optical lattice potential ltl 17 1. 
The kinetic energy and angular momentum per parti- 
cle are i?kin = (2/n)E ICC and L = 8H/ir 2 , respectively. 
The particle density p(r) tends to zero within the vortex 
cores, which have a radius on the order of i? CO rc ~ 1/k. 
Defining the healing length as £ = (87rap) -1 / 2 with 
mean density p, upon use of Eq. ([2]) one may write 
(^£) 2 = Erec/ 99 = E lec /V with mean potential well 
depth V. Thus, the condition that the core size exceeds 
the healing length £ < R colc is equivalent to E rcc < V. 

The stability of solutions ip of Eq. (Q} may be con- 
sidered via the grand canonical potential if^e] = 
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FIG. 1: (a) Flux density and (b) local particle density for 
sin(0) = 1. Black (white) indicates low (high) particle density. 
The dashed white square defines the position and size of the 
plaquettes used to build up the lattice. 



Jd 3 r[r £ f^ + (V tlap -^ e \ 2 +9\A\y2} for = ^ + e X 
with the kinetic energy operator T, e £ R, and an 
arbitrary normalized wavefunction \. Use of Eq. ^ 
yields f K[^] e=0 = and ^K[^ £ ] E=Q = f d 3 r[2 X *(f - 
E rcc )x ~ Hrap (e~ lS X + e lS X*) 2 ]- This integral is evalu- 
ated on each plaquette by expanding x in a Fourier series 
up to second order (higher order terms do not contribute 



because Vt iap (x,y) is of second order). Hence, by solv- 
ing an eigenvalue problem with finite dimensions we find 
that tjj is metastable in the sense that it establishes a true 
local energy minimum if V > 3I? rec . 

The experimental realization of I4rap as required by 
Eq. {2J may be achieved by implementing a light shift 
potential obtained from the superposition of two opti- 
cal standing waves arranged parallel to the x- and y- 
axes, with linear polarization parallel to the z-axis 18]. 
If g > 0, the trap potential satisfies Vt rap < and 
may thus be confining, which corresponds to adjust- 
ing a negative detuning of the optical standing waves 
with respect to the relevant atomic resonance. Available 
laser technology allows for lattice constants A/2 = n/k 
in a wide range between a few hundred nm and sev- 
eral /im. The d x 2_ y 2-wave solution ip does not rep- 
resent the ground state and thus, in addition to the 
lattice potential, an experimental procedure is required 
to excite it. In the following, a method based upon 
stimulated Raman scattering is discussed for the case 
of polarizable bosons in a two-dimensional square op- 
tical lattice. We extend an approach developed for 
vortex excitation such that we may apply angu- 

lar momentum quanta h with alternating signs to the 
A/2-sized square plaquettes. Let us consider a linearly 
polarized bichromatic light-field E(r, t) — z[E\(r,t) + 
E 2 {r,t)] defined by E x (r,t) = A x e iut \cf>(x, y)\e iS ^ and 
E 2 (r,t) = A 2 e i ^ +n)t \(t>{x,y)\e- tS ^^ with real posi- 
tive amplitude coefficients A\, A 2 , frequencies w, cj + 0, 
and \(j)(x,y)\ and S(x,y) as in Eq. ([3]). Here, |0| <C cj 
is assumed such that both fields, despite their differ- 
ent frequencies, in good approximation share the same 
cikonal S(x,y). The total intensity of the light field 
I(x,y,t) = E(r,f)E*(r,i) = I L (x,y) + I R (x,y,t) (aver- 
aged over the optical period 2ir/uj) consists of a station- 
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(a) 



FIG. 2: In (a) the square lattice intensity Il(x,ij) is plotted. 
White (black) regions indicate antinodes (nodes), which for 
negative detuning correspond to potential minima. In (b)-(j) 
the micro-rotor intensity In(x,y,t) is shown (same area as 
in (a)) for Sit = mr/8 with n G {0, 1, 8}. The grey scale 
indicates values of Ia(x,y,t)/AA\A2 between —1 (black) to 
+1 (white). The dashed white squares indicate the position 
and size of the plaquettes. 
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FIG. 3: (a) Optical set-up consisting of two nested Michelson- 
interferometers. PZT = piezo-electric transducer, M = mir- 
ror, BS = beam splitter, AOM = acousto-optic frequency 
shifter, (b) The lattice is decomposed into plaquettes (grey 
rectangles) with potential minimas in the corners (filled disks) 
and maximas in the centers (white disks) translated along the 
basis vectors indicated by the black arrows. 



ary term Ii,(x,y) and a time-dependent term In(x,y,t) 
with 



I L (x,y) = (A 1+ A 2 ) 2 \<f>(x,y)\' 



I R {x,y,t) = 2A 1 A 2 \<f>(x,y)\ 2 cos{2S(x,y)-nt). 



(5) 



Both terms give rise to light shift potentials according 
to Vl,r — — $l(ct)lL,R with atomic polarizability a [ioj ]. 
The amplitudes A\ and A 2 can be chosen to control the 
ratio Ir/Il- We restrict ourselves to sin(0) = 1 in the 
following. In this case, the stationary term provides a 
square lattice Ji(a;,y) = {A\ +^42) 2 [sin 2 (fcx) + sin 2 (A;?/)], 
while the time-dependent term Iji{x,y,i) acts as a col- 
lection of microscopic rotors, which apply angular mo- 
mentum with alternating sign to the A/2 x A/2-sized pla- 
quettes of the square lattice lL(x,y). As illustrated in 
Fig. (2), in the vicinity of each node of Il(x, y) the term 
Ifi(x,y,t) provides a rotating quadrupole potential with 
alternating sense of rotation for adjacent lattice sites. In 
Fig. 2(b)-(j) the micro-rotor intensity Ir(x, y, t) is shown 
for f2i = rar/8, with n £ {0, 1, ...,8}, thus illustrating a 
fit = 7r/2 clockwise rotation of the quadrupole around 
the central node of Ii,{x,y) in Fig. 2(a). Hence, the an- 
gular frequency of the micro-rotors is (l = Q/2. 

Temporary application of the micro-rotor potential 
Vft(r, t) may be used to resonantly excite the d x i_ y i- 
state tp by adjusting Ml equal to the energy difference 
between tjj and the ground state of the lattice 2l| . Let 



us briefly estimate the resonance condition. The angular 
momentum applied to each plaquette is approximately 
given by mSlr 2 , where r = A/4 is the distance from the 
centre to the edge of the plaquette. Excitation of vor- 
tices requires an angular momentum of h per plaquette, 
i.e., mtlr 2 h and thus hSl « (8/-7T 2 ) 2E rcc . For ru- 
bidium atoms and a convenient lattice wavelength (A = 
1030 nm) H/2tt — 3.5 kHz. Since Ml turns out to be com- 
parable to the collisional energy gp, we cannot directly 
apply the theory of Ref. 2l| to estimate the excitation 



efficiency. A detailed calculation is involved and beyond 
the scope of this article. Nevertheless, the appropriate 
transfer of angular momentum by Vr(t, t) together with 
the metastable character of i/j yield strong evidence that 
if) can in fact be efficiently excited by means of Vn(r, t). 

Experimentally, the generation of the bichromatic light 
field of Eq. ([5]) is straightforward using the optical set- 
up illustrated in Fig. 3(a), thus extending a method 
proven practicable in previous experiments 22]. The 
two components Ei(r,t) and E 2 (r, t) are produced in 
two nested Michelson-interferometers, each with its two 
branches folded under 90° angle. Two laser beams with 
adjustable frequency difference and linear polarization 
perpendicular to the drawing plane in Fig. 3(a) are used 
to couple both interferometers. Each interferometer com- 
prises a piezo-electrically driven mirror, which permits 
to adjust the appropriate relative phase of the fields in 
each branch. In order to produce the required waves 
(j){x,y) = e 471 "/ 4 sin(fca;) + e~ 47r / 4 sin(fcy) and 4>*(x,y) = 
e -iTr/4 s j n (/j2;) _|_ e W4 S i n (/ C y) j the optical path length dif- 
ferences are set to +A/4 and —A/4. 

Finally, we wish to point out that the micro-rotor term 
Iu(x,y,t) gives rise to driven ring-like tunneling in the 
Bose-Hubbard model. Expanding the boson field oper- 
ator in terms of the Wannier-function W(r) of the first 
band yields 



H R = J dr* + (r,t)I/ fl (r,t)*(r,t) 



Kb m J drW*(r-r n )V R (r,t)W( 



(6) 



where Vn{r,t) = V(r)e int + c.c. is the micro-rotor po- 
tential [dr = dxdy) and n, m are summed over all lattice 
sites. This yields H R = J2 n m b+b rn K nm e^ nt + H.c. 
with K nm = I drW*(r - r n )V(r)W(r - r rn ). Keep- 
ing only diagonal and nearest neighbour terms and as- 
suming that W*(r — r n )W(r — r m ) ^ mainly around 



4 



(r„ + r m )/2, one may write if nm 
with K Ilm = / drW*(r - r n )W{r 



m 

V((r n + r m )/2) 
satisfying n nm = 



K mn = K G R for n ^ m and k„„ = 1. The lattice can be 
composed by translations of A/2 x A/2-sized square pla- 
quettes (grey rectangles in Fig. 3(b)) along the Bravais- 
lattice TZ = {i?„ M = vR+ + fiR-\v,fi € Z} with basis 
vectors R± = (A/2) (a; ± y), respectively (black arrows 
in Fig. 3(b)). Assume that the four potential minima at 
each corner of some plaquette (filled disks in Fig. 3(b)) 
are indicated by n = 1, 2, 3, 4 starting from the lower left 
corner in clockwise order. Using V(r) = (Vi/2) <p 2 (x, y) 
with a real constant V\ , the coupling constants are eval- 
uated as K\2 — K 3 4 — —K23 — —K41 — iViK/2 and 
Kn = K33 = —K22 = = V\. Summing over all 

plaquettes finally yields 

H R = ViKsin{m)J2( b i -bi)(b3-h) + H.c. 



■R 



+ 2Vicos(m) (N+ -JV_) , 



(7) 



3 &3 and 



where the operators N + = (1/2) Y^n + & 
iV_ = (1/2) 6^6 2 + 64 &4 denote the total particle 
numbers in the two sublattices indicated by (+) and (— ) 
signs in Fig. 3(b). 

The Hamiltonian Hr has to be added to the conven- 
tional Bose-Hubbard Hamiltonian Hbh describing the 
optical lattice in absence of Vr(j", t) according to Ref. 0. 
The first term comprises next neighbour hopping terms 
oscillating in phase for opposite edges of a plaquette and 
with opposite phase for adjacent edges. The second term 
provides potentials oscillating in phase for lattice sites 
on the same diagonal and with opposite phase for ad- 
jacent sites. Hopping and potential terms exhibit 90° 
phase lag, thus giving rise to a rotational direction of 
Hr, which suggests the interpretation of a driven ring- 
like tunneling process. Note that Hr is quadratic with 
respect to the particle operators, in contrast to the quar- 
tic ring exchange interactions considered in Refs. Hill. 
The ground state of i?BH exhibits a well known phase di- 
agram possessing a Mott insulator and a superfluid phase 
8]. The implementation of Vp(r,t) allows to experimen- 
tally study the phase diagram of Hbh + Hr. Whether it 
provides additional quantum phases (e.g. a normal Bose 
liquid, which is compressible but not superfluid) is an 
exciting question to be addressed. Via temporary appli- 
cation of Hr one could explore the possible existence of 
an excited state of Hbk corresponding to the mean field 
solution given by Eq. |2]). 

In summary, we have discussed the existence of a 
d x 2 _ y i -density wave in a two-dimensional bosonic opti- 
cal lattice with square lattice geometry. We have pro- 
posed a straightforward experimental scenario to excite 
this state, which in the Bose-Hubbard picture amounts to 
the insertion of anisotropic time-dependent hopping and 
potential terms. This opens up a new axis of phase space, 



possibly giving rise to new exciting quantum phases. 
Our considerations are not constrained to square lattice 
potentials. For example, the superposition of six light 
beams travelling in the xy-plane (one every 60° with 
alternating frequencies) yields a similar scenario of al- 
ternating rotational flux, however with hexagonal lattice 
symmetry. In this article, we have restricted ourselves 
to bosonic atoms, however, similar scenarios involving 
fermions, mixed systems or molecules appear possible. 
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